Generalized Green-Kubo formula and fluctuation theorem for quantum current 
system around a nonequilibrium steady state 
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Fluctuation theorem is derived for a quantum current system around a nonequilibrium steady 
state. It is demonstrated that the fluctuation theorem can be a part of the generalized Green-Kubo 
formula or a nonlinear response theory of an external field or a change of the chemical potential 
difference. 
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Green-Kubo formula[T| is one of the most funda- 
mental relationships in nonequilibrium statistical me- 
chanics. The original derivation is restricted to the 
linearly nonequilibrium case, but there exist various 
generalizations. Recently, we have recognized that 
the fluctuation theorem is a nonlinear generalization of 
Green-Kubo formula, and is closely related to some iden- 
tities such as Jarzynski equality. p7H13| The fluctuation 
theorem can be applied to various fields and is impor- 
tant in particular for mesoscopic physics. 

So far, the fluctuation theorem may be regarded as a 
nonlinear response theory around an equilibrium state. 
Indeed, one of the key steps is to derive the fluctuation 
theorem is to assume a local detailed balance condition. 
Thus, nonlinear response theory around a nonequilibrium 
steady state is outside the existing fluctuation theorem. 

The purpose of this paper is to generalize the fluctua- 
tion theorem around a nonequilibruim steady state. The 
formulation is a quantum extension of works by Chong 
et al. 14 1 and by Evans and Morriss@ as well as a pre- 
liminary formulation by the present author. (l5l] 

Let us consider a system coupled with two heat baths 
in which one of the chemical potentials is higher than an- 
other. The total Hamilitonian H to t(t) consists of three 
parts; the system Hamiltonian ifg(i), the bath Hamilto- 
nian Hb and the interaction Hamiltonian Hi nt (t) as 



fltot(t) = H s (t) + H B + gH int (t), 



(1) 



where g is the coupling constant. 

The setup we consider in this paper is as follows: At 
first, we prepare an equilibrium system as the initial con- 
dition in which there is no coupling between the system 
and the bath, i.e. H- mt (t) = at t — > -co. Then, we wait 
for the system to reach a nonequilibrium steady state af- 
ter an adiabatically connecting between the system and 
the bath until t = 0. At t — 0, we may disturb the system 
and investigate a response to this disturbance. 

For later discussion, let us consider a nonlinear re- 
sponse around a nonequilibrium steady state obtained 
from a coupling between the system and the bath. Equa- 
tion (JT|) is the most general form of the Hamiltoninan, but 
it is difficult to extract some useful information based on 
this treatment. Indeed, a density matrix for a nonequilib- 
rium steady state is necessary for the description of heat 



conduction and the electric conduction under a chemi- 
cal potential difference, but it seems to be impossible to 
obtain the steady density matrix for Eq. (JJ . It is remark- 
able, howeverjthat Fujii[f| derived a steady McLennan- 
Zubarev form[3|,[l6[ of the density matrix under a special 
setup within the framework of Keldysh formalism. [17H19l | 
Here, we outline the derivation of the McLennan- 
Zubarev density matrix within the framework of Keldysh 
formalism. Here, let us consider the case that the reser- 
voirs is adiabatically connected with the system as 

Hint, 



H{ ot (t) = H + ge 



(2) 



where Hq = Hg + Hb under no time dependence in Hq (t) 
and Hint(t), H* ot (t) — > H tot (t) in the limit e — > 0, and 
the coupling between the system H(t) and the bath H-q 
is adiabatically zero in \t\ — > oc. Note that we always 
take the limit e — >■ for later discussion, even when we 
do not write its limit explicitly. 

Here, the density matrix at t = is given by 



p e = S e (0,-oo)p°S e (-oo,0), 



(3) 



where S e (t,t ) = 7Vexp[-J- 5 dse~ e ^ H int (s)] for t > 

t Q and S e (t,t ) = 7%cxp[--|5 dse~ e ^ H int (s)] for t < 
to with the introduction of the right ordered operator 
and the left ordered operator T<_. Here, the initial 
condition of the total system in the limit t —> -co is 
assumed to be 



-l3(H -e$AN/2] 



/Z t , 



(4) 



where e is the charge of an electron, $ is the voltage dif- 
ference between two reservoirs, and A7V is the number 
difference of electrons between two reservoirs. To repro- 
duce this initial density matrix (0]) we have assumed that 
the average (/^l + Mr)/2 of the left chemical potential /xl 
and the right chemical potential /xr is zero, and 4> satis- 
fies e$ = /il — MR- Note that the expression ^ is only 
valid for the case of [Hq^d ] = 0. 

As was shown in Ref. p , it is notable that any observ- 
able A commutable with H at t = satisfies 

r° 

A e = S e (0,-oo)AS e (-oo,0) = A+ / dte- e ^jfi(t) 

J — OO 

= e f dte-^A H (t), (5) 
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where we have introduced jfi(t) = -{d / dt)A H (t). Note 
that A e in Eq. ([5]) is equivalent to the invariant measure 
introduced by Zubarev.Q 

It is straightforward to obtain the density matrix at 
t = as@ 



p e (0) = exp \-p 



if £ (0)-^A7V-(0) 



/Z tot (6) 



in the limit e — > 0, where 



i^(0) = H + / ^e- £ l*l J e , H (t), 

■/ — OO 

/.O 

AiV e (0) = AN+ dte-^J cM (t) (7) 



with the energy current J* H (t) = —(d/dt)Ho,H{t) and 
the mass current J* H (t) = — (d/dt)ANH{t). The density 
matrix (JSJ is equivalent to the classical nonequilibrium 
density matrix by McLennan- ZubarevJH For later 
convenience, we rewrite Eq.fJH]) as 



P e (0) 



exp{-^ e (0)} 

zW) ' 



(8) 



where H e (0) = H^O) - e($/2)AAT £ (0) and Z e (f3) = 
tre -/m E (o)_ Thig ig the 

consequence of Keldysh Green 
function formalism for a nonequilibrium steady state. Q 
It is straightforward to introduce p e (t) at t satisfying 
\t\ < 1/e as the generalization of Eq.(|6]). It should be 
noted that the density matrices p\ ot (t) is in a steady state 
in time \t\ <C 1/e. 

We should note that von-Ncumann equation for the 
total system 



dt 



p\t) = -iC\ otfi {t)f{t) = ±{HUt),p e (t)} (9) 

in 



for the density matrix p\ ot {t) satisfies 
dp e (t) 



dt 



iC\ ot Jt)?{t)=Q 



(10) 



for \t\ <C 1/e, where the subscript for the Liouville oper- 
ator C%{t) is introduced to represent a reference state for 
the response theory for later discussion. Because of the 
relation [p £ (0), # e (0)] = [p e (0), AiV e (0)] = the density 
matrix p e (0) also satisfies 



iCt ) (t)p*{t)=Q 



(11) 



for \t\ <C 1/e, where the Liouville operator C^{t) can be 
formally written as 



iC e (t) = -[W(t), } 



(12) 



for \t\ <C 1/e. It should be noted that the Liouville oper- 
ator £qW i s unitary. 

Now, let us discuss how we can construct a nonlin- 
ear response theory around the nonequilibrium steady 



state Eq.©. To avoid heavy notations, we drop the adi- 
abatic factor e of any variables for later discussion un- 
der the limit e — > 0. Therefore, though our argument 
is only valid < t < 1/e, the formulation can be used 
for any finite t in the limit e — > 0. We should note that 
iCo(t) = lim e _>o *£q(0) for \t\ <C 1/e. We consider re- 
sponse theories after the following two cases for t > 0: 

= $ + A$(i), (case A) (13) 
H (t) = H (0) + F cx (t)B (caseB) (14) 

where B is the conjugate field to the external force F cx (t). 

Let us introduce the difference of the Liouville opera- 
tors iA£(t) = iC(t) — i£o(0) from the base state as 



iA£ (A )( , ) = _m(!) [AiV! 



for case A and 



iAC^\i) 



iF ex {t) 



[B, 



(15) 



(16) 



for case B, where iC(t) is the Liouville operator for t > 0, 
and the upperscripts (A) and (B) are introduced corre- 
sponding to the cases A and B, respectively. (To repre- 
sent both cases we do not put any upperscript) . By using 
iA£(t), the time evolution of the density matrix of the 
system at t > can be described by 



-p(t) = -iA£(t)p(t), 
and its formal solution is given by 

p(t) = U<-(t,0)p(0) 

where 

£/<_(*, 0) = r < _[e -i J'o 



(17) 



(18) 



(19) 



These results are based on Schrodingcr picture. 

We can also introduce the adjoint dynamics for observ- 
able Ajj(i) based on the Heisenberg picture as 



j t A H {t) = U^(-oo,t)i£{t)A, 



(20) 



where the subscript H represents an observable in Heisen- 
berg picture, and we have introduced 



U^{t 0l t)=T^e lf ^ dsC{s) . 



(21) 



(22) 



It is possible to prove the relation^ 

U^(t U t 2 ) = U^{t U T)U^(T,t 2 ) 

and 

U^(t ,t) = U^(t,t )- 1 . (23) 
The formal solution of Eq. ([20| is given by 



A H (t) = U->(p,t)A H (P) = U^(p,t)A H (P)U^(t,0), 

(24) 
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where A H (0) = U%(-oo,0)A = T^e 1 S -™ dsCo[s) A. The 
last expression can be proved with the aid of Baker- 
Campbcll-Hausdorff formula. Here, E/_>.(0,i) in Eq. (|2"Tj) 
can be rewritten as 

(0,t) = T_>eJ? dstAC( -^ , (25) 

because ^4_y(0) also satisfies the steady condition 

E7^(0,t)A(0)=i4(0) (26) 

with A(0) = A H (0) and ?7°.(0,t) = T.+eti dstC °( s \ Note 
that there are two relations 

B H (t) = B, and AN H (t) = AN. (27) 

Therefore, the expectation value of the observable A 
based on the steady density matrix p(0) defined in Eq. 
© is represented by 

(A H (t)) = tT S {p(0)A H (t)} =tr{p(p)U->(O t t)A H (P)}. 

(28) 

The Dyson equation between U-y(0,t) and U^.(0,t) can 
be written as 



,(0,t) = C7* (0,t) + f dsU%(0,s)iA£(s)U^(s,t). 
Jo 

(29) 

From the steady condition (|26|) we obtain 
(5A H (t)) 

dTtT{p(0)U°,{0,T)iAC(T)U^(T,t)A H (0)} 

drtr { [Ul (r, 0)p(0)] i A£(r) 17_> (r, i) A ff (0) } 

drtr {p(0)iAC(T)U^(T,t)A H (Q)} , (30) 

where ((L4jj(t)) = (Au(t)) — (Ah(0)), and we have also 
used the equivalency between Schrodinger picture and 
Heisenberg picture in the second line and the steady 
condition U^_(t, 0)p{0) = p(0) for the final expression. 
Equation (|30|) can be rewritten as 



(SA H (t))= / dTtr{[iAC(T)p(0)]U->(T,t)A H (0)} 
Jo 

(31) 

in terms of the equivalency of two pictures. 

Because of the relation ([6]) and (|27p , Eq. ([T5"T) is rewrit- 
ten as 

iAU(t)p(0) = - l -^§P-[AN H (t),p(0)} 
eA$(f) rf) 



where we have used Kubo's identity in the second line 
and we have introduced {t) = l/(ih)[AN H (*), H] 
and f H N \t;T) = e -«*(°W» (t)e- ffl (°W 6 . Similarly 
Eq. tfTTj]) can be rewritten as 



iA£ B (t)p(0) = F cx {t)p(0) / dXJ B (t;-ihX), (33) 



where J^'it) 



l/(ih){B H {t),H(0)] and J£\p,t) = 

e -iU(0)r/hj(B) ^ e -iH(0)r/h_ 

Because case A and case B are almost equivalent with 
each other because of the comparison between Eqs. (|3"2")) 
and (|3"3")l . it is sufficient to prove one of two cases. In this 
paper, let us analyze the case B. Substituting Eqs. (fTTj) 
and (|33[) into Eq. (|3"Tj) we obtain 



(6A H (t)) 



[ drF cx {r) t ' d\ 
Jo Jo 

(J ( *\T;-ihX)U^(T,t)A H (0)). (34) 



With the aid of C/_>(r,t) = f7_,(0, r)- 1 J7^(0, t) 
[/<_(-?-, 0){7_>(0,i) Eq.([53)l is further rewritten as 



(6A H (t)) 



dTF cx (T)(J i I f > (T;-ihX)U^(s,0)A H (t)} 



- [ drF^r) [ 
Jo Jo 



dAtr{C7_,(0,r) 



p(0)jg\<r,-ih\) A H (t)} 



(35) 



With the aid of Eq. (f27j) we have the invariant relation 
J y H ' (t) = J H . Then, we may have the relation 

U^(p t tM0)J^\t;-ihX)] 

= U^(0,t)p(0) ■ U^{0,t)J ( I f\t;-ihX) 



= U^(0,t)p(0) ■ J { H ] {t;-ihX). 
Here, from Eqs. ([T9l and (I3"3"|) we can write 

U-y{0,t)p(0) =p{Q)T^e-ti drn{ - T \ 



(36) 



(37) 



where f2(r) = — F ck (t) Jq dXJ^ (t; —ihX). Substituting 
Eq.(j37|) into §3S§ we obtain 



p(0) / dXe x ^J { H N \t)e- xn ^ 



{5A H (t)) = / dsF 0X (s) / dX^ BA (s), (38) 
Jo Jo 

where = {T^e-fo dTU ( T ) jjf \t; -iHX)A H (t)) is 

the nonlinear response function. Equation (|38l) is 
the quantum version of the generalized Green-Kubo 
formula. If we replace the observable A by the current 
J, the expression ([55} can be rewritten as 



e -^-p(0)J o dXJ^it-iKX), (32) 



(<5J ff (i)) 



/" dsF ex (s) [" d\$ B j((a), 
o Jo 



(39) 
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where we have introduced {SJnit)) = (JH(t)} — (J(fy)- 
These expressions are the quantum version of the non- 
linear response theory obtained by Evans and Morriss. Q 
Needless to say, it is reduced to the conventional linear 
response formula[l| if we set Q(t) = 1. It should be noted 
that the nonlinear response theory in case A can be ob- 
tained when we replace J^\t; —ihX) by J^^r; —ihX). 

The fluctuation theorem is easily obtained from the 
normalization condition trp(-t) — 1 as 

(T^e- X> drn(T) ) = 1, (40) 

where we have used p{— t) = U-^(0, t)p(0). Eq uation (|4TJ)) 
is the integral fluctuation theorem in Refs.fliTflij. From 
Jensen's inequality, (j40|) leads to J Q dr(17(T)) > 0, which 
means that gentropy production" is positive through the 
process. 

It should be noted that our nonlinear response theory 
of fluctuation theorem has the same form as that around 
an equilibrium state thanks to the special choice of the 
nonequilibrium steady state Eq. (|6]) . Of course, we cannot 
obtain the explicit form of any expectation value within 
the framework presented here, because the steady density 
matrix Eq.® contains the history of the current. The 
explicit calculation to demonstrate the usefulness of this 
framework will be reported elsewhere. 



In conclusion, we have derived a quantum nonlinear re- 
sponse theory (|3"5)l around a nonequilibrium steady state 
([6]) for a dissipative quantum system. We have also de- 
rived the quantum fluctuation theorem (j40[) for dissipa- 
tive quantum systems. This formulation is directly ap- 
plicable to the description of the dynamics of a subsys- 
tem charcterized by the system Hamiltoninan Hs if the 
coupling between the system and the bath is sufficiently 
weak. Actually, Evans and Morriss Q derived the gener- 
alized Green-Kubo formula for classical dissipative sys- 
tem under the assist of the adiabatic incompressibility of 
the phase space. The formulation should be extended to 
the case for dissipative subsystems with finite coupling 
between the subsystem and the bath. Although the re- 
sults presented here is rather formal, the explicit calcula- 
tion following this formulation will be reported elsewhere. 
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